TTP04-20 
SI-HEP-2004-09 



Modeling the Pion and Kaon Form Factors in the 

Timelike Region 



Christine Bruch a \ Alexander Khodjamirian h \ 
Johann H. Kiihn a ) 

a ) Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, 

D-76128 Karlsruhe, Germany 

b ) Theoretische Physik 1, Fachbereich Physik, Universitat Siegen, 

D-57068 Siegen, Germany 



Abstract 

New, accurate measurements of the pion and kaon electromagnetic form factors 
are expected in the near future from experiments at electron-positron colliders, using 
the radiative return method. We construct a model for the timelike pion electromag- 
netic form factor, that is valid also at momentum transfers far above the p resonance. 
The ansatz is based on vector dominance and includes a pattern of radial excitations 
expected from dual resonance models. The form factor is fitted to the existing data 
in the timelike region, continued to the spacelike region and compared with the mea- 
surements there and with the QCD predictions. Furthermore, the model is extended 
to the kaon electromagnetic form factor. Using isospin and SU(3)-flavour symmetry 
relations we extract the isospin-one contribution and predict the kaon weak form 
factor accessible in semileptonic r decays. 



1 Introduction 



The pion electromagnetic (e.m.) form factor F n (s), one of the traditional study objects 
in hadron physics, nowadays plays an essential role for the precise determination of elec- 
troweak observables. An accurate knowledge of F n (s) at timelike momentum transfers 
s > 4m^ is needed to calculate the hadronic loop contribution to the muon anomalous 
magnetic moment and to the running of the e.m. coupling (see |2| for the current 
status). 

New accurate data on the pion form factor in the timelike region have recently been 
obtained by the CMD-2 collaboration [3| measuring the e + e~ — > 7r + 7r~ cross section at 
yfs = 0.61 -r 0.96 GeV (for an update of these data see 0j). In this region, it is successfully 
described (fitted) using models based on p-meson dominance, with a small but clearly 
visible cu-meson admixture. Above 1 GeV data on e + e~ — ► n + 7i~ jSUHlIZI exist but are not 
that accurate. Employing isospin symmetry one also gains independent information from 
the measurements of r — ► 7r _ 7r z/ T at s < m 2 (for details see, e.g. PQ). There are other 
interesting form factors closely related to F n : the charged (neutral) kaon e.m. form factors 
measured in e + e~ — > K + K~ (e + e~ — > K°K°), as well as the weak transition form factor 
accessible in r — > K~K°v T . 

In the near future the experimental knowledge on F^^is) will be substantially im- 
proved, due to new data to be obtained using the radiative return method jj.8.. The first 
measurements of F v with this technique in the p region have already been performed by 
the KLOE Collaboration and the agreement with the CMD-2 data is encouraging. At 
larger energies, up to 2.0-2.5 GeV, perhaps even 3 GeV, accurate measurements of F W)K 
are anticipated from the BABAR experiment (for preliminary results see [TO]). The high 
rates expected at s > m 2 p demand phenomenological models more elaborated than the 
simple p- (p,u,(j)-) dominance models for F^ {Fx). The main purpose of this paper is to 
construct an ansatz for the pion form factor that is valid in the region below and far above 
p-resonance and to extend this model to the kaon form factor. 

The model is constructed to obey the constraints from analyticity and isospin-symmetry 
and to incorporate the proper behaviour at high energies, consistent with perturbative 
QCD, and the correct normalization at s — 0. Furthermore it is based on plausible as- 
sumptions derived from the quark model, moderate SU(3)-breaking, vector dominance and 
a pattern of radial excitations expected from dual resonance models. It has enough flexi- 
bility to accommodate the characteristic interference pattern of the cross section and, once 
sufficiently precise data are available at higher energies, may be used to fix the parameters 
of the higher excitations. 

Above the p-resonance the excited p'(1450) and p"(1700) are expected to play an impor- 
tant role. Already now these two states are indispensable, if one wants to accommodate 
the measured parameters of the p-resonance with the correct normalization of F w (s) at 
s = 0, as demonstrated by earlier analyses and fits within the p region (see, e.g. fH\ 1X2]). 
In general, to include all possible intermediate hadronic states in the 7* — > 7r + 7r~ transi- 
tion amplitude, one has to take into account an infinite series of radially excited p's. In 
addition, there are multihadron intermediate states with J p = and 1 = 1 (2ir , An , K K 
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etc.). Hence, the pion form factor at timelike s > m 2 p is a complicated object determined 
by a large or even infinite amount of hadronic parameters not accessible at present in a 
rigorous theoretical framework. 

On the other hand, at sufficiently large s one expects F„-(s) ~ a s (s)/s, as predicted 
from perturbative QCD ^3] for the pion form factor in the spacelike region at s —>■ — oo, 
analytically continued to s —>■ +00. In other words, the overlap of many intermediate 
hadronic states has to build up a smooth, power-behaved function. One might consider 
using this QCD prediction in the region of present interest, that is, at intermediate timelike 
s. However, data on the form factor in the spacelike region, s < 0, indicate that the 
onset of this asymptotic behaviour is far from the "few GeV 2 " region. Preasymptotic 
contributions ~ l/s n with n > 1, stemming from the end-point, soft mechanism are 
essential for s ~ 1 — 10 GeV 2 . Approximate methods valid at intermediate spacelike 
momenta, for example QCD sum rules jT^l EH UZj , allow to calculate F n (s) including soft 
effects. However, a straightforward analytic continuation of F n (s < 0) to large s > is 
difficult. The timelike form factor will suffer from uncertainties in the analytic continuation 
of soft parts, Sudakov logs and of a s (s) (for a discussion see [TBI EH])- Hence, QCD 
calculations cannot be directly used in the large s > region, e.g., for estimating the 
"tail" of higher resonances in the form factor. In this paper we therefore prefer to adopt a 
model for the pion form-factor formulated entirely in terms of hadronic degrees of freedom. 
Importantly, the form factor, as obtained from the model, fitted in the time-like region 
and extrapolated to spacelike momenta, has a proper power-law behaviour which can be 
compared with the data at s < and used to test various QCD-based predictions. 

Models of hadronic amplitudes, where an infinite series of resonances at s > is summed 
to yield a power-law behaviour at s < 0, are rooted in the Veneziano amplitude and dual 
resonance models formulated long before the advent of QCD. Importantly, the pattern of 
infinite zero-width resonances is predicted in the N c = 00 limit of QCD. Recently, the model 
for the pion form factor using the masses and coefficients chosen according to Veneziano 
amplitude was considered in [20]. Earlier, similar analyses of the pion form factor can 
be found in |2T] . Models of dual- resonance type with infinite number of resonances are 
widely used also for other hadronic problems, some recent works can be found in [22] . 
We will use the dual-QCDAr c=00 model j20] as a starting point, modifying it for the first 
few p resonances, by keeping their parameters (masses, widths and coefficients) free and 
fitting them to experiment. In this way, the complicated effects of p resonances coupled to 
multihadron (27r,47r etc.) states are implicitly taken into account. It is remarkable that 
the gross features of the model are well reproduced with the (fitted) resonance parameters. 
Since in the dual-resonance amplitude the coefficients of higher resonance contributions 
decrease with the resonance number, the corresponding modifications for individual higher 
states are not important, and the "tail' of resonances is treated as in J20]. The model for 
the pion form factor is also analytically continued to the spacelike region and compared 
with the data there and with the QCD predictions on F n (s) at large spacelike s. 

Furthermore, we extend the model to the kaon form factor, employing an SU(3)- 
generalization of the pion amplitude. Fitting the charged and neutral kaon form factors 
to the data, and using flavour symmetries, we predict the weak kaon form factor relevant 
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for r semileptonic decays. Let us also mention at this point that the decomposition of the 
form factors into their isospin zero and one components respectively, is also of relevance for 
a model-independent evaluation of 7-Z-mixing j2Hj where the two amplitudes contribute 
with a different relative weight. 

The plan of the paper is as follows. In Section 2 we summarize the phenomenology 
of F n (s) recalling the derivation of the p-meson contribution, whereas in Section 3 the 
contributions of excited p resonances are discussed. The model for the pion form factor is 
introduced in Section 4 and its parameters are fitted to the data. In Section 5 we proceed 
to the kaon e.m. form factors and Section 6 is devoted to the weak kaon form factor in r 
decays. Section 7 contains our summary and conclusions. 



2 The p-meson contribution to the pion form factor 

The pion e.m. form factor is defined in the standard way, 

(vr + (pi)7r-(p 2 ) I f™ I 0) = (p! -p2)Ma) (1) 
The quark e.m. current = Yl e q ( ll^ ( l can be decomposed into isospin one and zero 

q=u,d,s 

components respectively. At the quark level this corresponds to 

•em = J_ .3 , JL -j=o _ I as (2) 
\pl 3^/2 3 

with 

jl = («7a» w - d^d)/V2 , j J =° = (ujf.u + V2 , = s7 M s . (3) 

The isotriplet partners of the current form the charged weak current 

= (jj + ijJ)/V2 = ti7 M d. (4) 

In the isospin symmetry limit, the 1 = and s-quark components of the current do not 
contribute to F n . In Eq. (jTJ), s = {pi + P2) 2 is the timelike momentum transfer squared, 
s > 4m^. The form factor F^(s), being analytically continued to the spacelike region s < 0, 
corresponds to the hadronic matrix element (7i + (pi) \ j^™ 1 \ 7r + (— p 2 )) related to Eq. ((IJ by 
crossing-symmetry. 

There are very few model-independent relations determining or constraining the pion 
form factor. One of them is the normalization condition for the pion electric charge, 

^(0) = 1 • (5) 

An important role is played by the dispersion relation, 

FM= iL^m. ( 6) 

n J s' — s — ie 

4ml 
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The asymptotic behaviour expected from perturbative QCD [T3], 



a 



lim F^s)-—, (7) 



s 



allows for an unsubtracted dispersion relation. The application of Eq. (jHJ) is based on an 
independent equation for the imaginary part of the form factor derived from the unitarity 
condition, 

2ImF jr ( a )(p 1 -ft), = £ J dT h (7r + ( Pl )7r-(p 2 ) \ h)(h | j™\Q)* , (8) 

h 

where all possible hadronic states h with J PC (I ) = 1 (1 + ) are inserted. Each term in 
the sum in Eq. (JHJ) includes the integration over the phase space and the summation over 
polarizations of the intermediate state h. Since isospin symmetry is not exact, there could 
also be a small admixture of the isospin-zero J PC = 1 states, e.g. the to and its radial 
excitations. 

Experimental data reveal that at low s < 1 GeV 2 the most important contribution to 
Eq. (jHJ) stems from p meson. The p-meson decay constant: 

<P° I 37 I 0) = , (9) 

and the strong pTin coupling: 

(vr + (pi)7r-(p 2 )|p ) = (pa - pi) a e^g^ , (10) 

(where e^ p ' is the p-meson polarization four-vector) determine the p-contribution to the 
imaginary part of the pion form factor in the narrow width approximation: 

ImF^\s) = ^g p7T7T n5(s-ml). (11) 

Substituting this into the dispersion relation (jHJ) gives: 

FM( 8 ) = ™f 9t ™ . , ■ (12) 
V l[m z p — s — te) 

The excited p', ... resonances have contributions of the same form, with the decay constants 
fpi^ and strong couplings (jy,™-, .... 

The pion form factor constructed by adding up the zero- width p, p', ... -resonances is an 
oversimplified ansatz which cannot be used at s > where the experimentally observable 
large widths of these resonances are important. The widths are generated by the contribu- 
tions of multi-hadron states to the imaginary part of F n , starting from the lowest two-pion 
state. The contribution of the latter to the unitarity relation 

2ImF^( S )(p 1 - p 2 ), = J dr 2 M - p^A^F^s) , (13) 



4 



involves the two-pion phase space : 

= J^2E[ (2^ (27r)454 ^ + ^ " ^ " ^ 

where (p'j + p' 2 ) 2 = s and A^s) is the amplitude of the strong pion-pion P-wave elastic 
scattering: 

A^(s) = (n + (p 1 )n-(p 2 )\n + (p' 1 )n-(p' 2 )) I=1J=1 . (14) 

In the low-energy region 4m 2 < s < 16m 2 , the two-pion state is the only contribution 
to the unitarity relation 1 . At larger s, vector mesons and various multihadron states 
contribute, making a model-independent use of Eq. (jHJ) impossible. 

A well-known and experimentally supported approach which we adopt here is "vector 
dominance" . In its simplest form it assumes that the p resonance saturates the pion form 
factor (Eq. (fT2|)). thus requiring j ~ p g p7r7r / ' \plm p = 1. One furthermore approximates the 
pion-pion scattering amplitude with an intermediate p exchange. Inserting the intermediate 
p propagator in Eq. (fPfj) and using the definition (fTUj) we obtain: 

A ( ^ (vr + ( Pl )7r-(p 2 )|p )(p |7r+(K)vr-(^)) 9%*(P2 - Pi) • (p' 2 

m p — s [nip — s) 

Substituting this ansatz together with Eq. (JT2j) in Eq. ([T3j). and integrating out the two- 
pion phase space with the help of 



J dT 27r (p' 2 -p' 1 ) a (p' 1 -p' 2 )n= (gan 



{Pl+P2) a {pi+P2)^\ [P{s)f , lg . 



where p(s) = |(s — 4m 2 ) 1 / 2 is the pion momentum in cm. frame of two pions, we finally 
transform Eq. (JTHJ) to 

ImF^(s) - mpfp f &MW ) ggZEZE ( 17 ) 
K [) ~V2(ml-s)\ 6nV~s ) m 2 p - s ' W 

This formula justifies using the more general expression 

F^)( s ) = m pfp A (2n] ("^ 9pn7T 



s 



V2{m 2 p -s) w m2- s ' 

which can be interpreted as a two-pion loop insertion in the p meson propagator. The 
amplitude A^^s) has a real and imaginary part. A natural approximation for ImA^ 
is the expression in curly brackets in Eq. ([170. At s — m 2 it is normalized to the p — > 2n 
width 

lmA^\m 2 p ) = m p T{p - 2vr) = [p{m 2 p )f , (19) 



For a review on the low-energy pion form factor see |24|. 
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To account for all possible amplitudes with two-pion insertions in the p-meson propa- 
gator, Eq. (fTHj) has to be added to Eq. (JT2J), together with subsequent diagrams with two, 
three, etc. two-pion loops. Summing up this geometrical series yields: 



rn 



2 



\f2mj m 2 p - s- ReA^is) - UmA^ (s 



for the part of F n which contains, in addition to the p-meson, the contributions of the 
2-pion intermediate state coupled to p. 

Several options for Eq. (j2*U|) are in usage. The simplest one is to neglect the real part of 
A^ 2n '(s) and to approximate the imaginary part by a constant, given by Eq. (|19|). This gives 
the usual Breit-Wigner (BW) formula for the p resonance with a constant width. A more 
refined version (used e.g. in |12j ) takes into account the s-dependence of Im^4^ 27r ^(s) in the 
form of the p-wave two-pion phase space (as indicated by Eq. (|I7j)) with the normalization 
from Eq. (JTHJ): 



ReA^(s) = 0, 1mA™ {s) = y/l^ ( T(p -> 2tt) = yGT^s) . 

s \p{m 2 )J 



(21) 



The function T p ^2n{s) naturally vanishes at s < 4m 2 , below the 2% threshold. 

The Gounaris-Sakurai (GS) approach [23 represents another option widely used. In 
this case one takes into account a nonvanishing real part of A^\s) calculated from the 
dispersion relation with two subtractions at s = 0: 

Using the expression for the imaginary part given in Eq. (|2*T|) and changing the integration 
variable from to = 4m^/sO one transforms the integral in Eq. 



s' z {s' — s — te) \ S[p{m p )\ 6 



4mJ 

1 1 

(23) 



I(s) = - I dv'v' A 



V 







v' — v — ie v' + v 



Calculating the principal value of the integral yields the real part: 

2 . „ 2 3, l + V 



Rel(s) = s ( - + 2v 2 - v 6 log — . (24) 

Using the above result in Eq. (|2*2*|) one obtains the real part of the amplitude A^ 2n \ Fur- 
thermore, following [25J the subtraction terms are fixed by the normalization conditions 
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for the mass and the width of the p: ReA i2w) (m 2 p ) = and £ReA {2n) (m 2 p ) = . At s = 
the form factor is normalized to unity. Hence 

w U \V2mJ m 2 -s + H{s)-i^T p ^{s)' 1 1 

where we use the same notation as in [T2] : 

H{s) = H{s) - H(m 2 p ) - (s - m 2 p )j-H{m 2 p ) , (26) 

so that 



H(s)= ^-&-^r ){-,-ml)v\og^. (27) 



'm 2 T(p -> 2n)\ f S 2 \ , 1 + v 
27r[p(m p )] J 

From experiment r(p — > 27r) ~ T tot (p), hence the couplings of p to other intermediate 
states can be safely be neglected. Therefore we replace in both versions of the B W formula 
the p — > 2tt width by the total width removing the superscript 2n at the form factor. 

Finally, the p contribution to the pion form factor introduced in Eq. (|T2*|) in the zero- 
width approximation and modified to include the width in Eq. ([20)1. can be rewritten in 
the following generic form: 

F^(s)=c p BW p (s) (28) 

where c p = F^ p \o) is the normalization coefficient. In the adopted approximation is 
determined by the product of p decay constants and pint coupling: 

c p = , (29) 

y 2.m p 

and BWp(s) is the BW formula normalized to unity at s = 0. For this formula two different 
versions will be used, one taken from 
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771 

BW p KS (s) = — p -=— - , (30) 

and the one from |2*H] 

m 2 + H(Q) 

w °« = ^F7T%fW)- (31) 

In both cases the effective s-dependent width is assumed to be 

r p (s) = T p ^ 27T (s) , (32) 
with r.h.s. defined in Eq. (J2T|) and T p (m 2 ) = T p ot . 
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3 Contributions of excited p states 



As already realized in earlier analyses of the pion form factor (e.g., in fUEj), the single 
p- meson approximation is not sufficient to fulfil the normalization condition (J5)l . Indeed, 
taking the measured values for T(p — > l + F) and T(p — > 2ir) from [2B| we obtain f p = 220 
MeV and, respectively g pn7T = 6.0, yielding c p ~ 1.2. One needs to include the contributions 
of excited p resonances to restore the correct normalization. Currently, two of them, 
p'(1450) and p"(1700), are well established experimentally [2E]. Adding the contributions 
of these two states in the form (J28|) to , one fits experimental data on e + e — > 2n in 
the p-region, practically up to s = 1 GeV, restoring the normalization condition F n (0) = 1. 
Both models (jHOj) and (jHljl work well. In addition, there is a small isospin-violating effect 
from u noticeable in the vicinity of p. In what follows, it will be taken into account as in 
[T2"] , by adding a. p — u mixing term to the p contribution: 

F ( P ) {s) ^ c j^M {1 + CwBWuj y (33) 

Note, however, that the dominant decays of excited p's are to final states other than 
2tt, hence the couplings of these resonances to various multiparticle intermediate states 
(Arc, KK etc.) become important at larger s, the region of our interest. In the previous 
section we have seen that the coupling of p to the 2tt state results in a geometrical series of 
27r-insertions into the p propagator yielding an imaginary part normalized to the p — > 2tt 
width in the formula for BW p . The analogous summation procedure can be repeated for 
each multihadron state coupled to a given excited p, say, to p'(1450). This leads to a 
formula for the p' contribution to F n similar to Eq. ()28|) . where p — > p' and the effective 
width in BW p i(s) is a sum over the effective widths for each channel 

T p ,{s) = T p ,^ 2 A s ) + r p' ->47r(s) + rV_* p27r (s) + T p ,^ KR (s) + ... (34) 

All we know about the functions on r.h.s. is their normalization at s = m 2 p , to the corre- 
sponding partial width of p', so that altogether Y p i{m 2 pl ) = r*°*. The s-dependence for each 
partial width has to be introduced in a model-dependent way, requiring detailed informa- 
tion on 47T (see e.g., j2Z]) and other hadronic final states in e + e~ . In particular, to account 
for a proper threshold behaviour one has to introduce phase space factors for each in 
Eq. ()34|) . different from the p-wave factor for T p (s). A complete kinematical and dynamical 
analysis of the partial widths for excited p resonances is beyond our task (some models 
can be found in [2EJ). We will continue using the same ansatz as for p, with the simple 
p-wave threshold factor, having in mind that there is still a room for improvement at this 
point. We have checked that small modifications of the effective width, e.g., replacing the 
effective threshold by Am^ have little influence on the form factor. 

Furthermore, the couplings of different vector resonances to one and the same multi- 
hadron state generate mixing between these resonances. To give an example of this effect 
return to the unitarity relation (J13)) and substitute on r.h.s. the p' resonance contribution 
to the form factor while keeping the intermediate p exchange for A nn . This term corre- 
sponds to a chain of transitions 7* — > p' — ► 2-n — ► p — > 2-n. This nondiagonal amplitude 
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is clearly not accounted for by the 2ir insertions to the individual p- and p'-propagators. 
One has to add to the form factor new terms with the products of two BW propagators, 
e.g. the p' contribution to the form factor will have the following schematic form: 

Fy\s) = c p ,BW p ,(s)(l + x p , p (s)BW p (s) + ...) (35) 

where Xpi p is the mixing amplitude which is s-dependent, in general, and ellipses indicate 
the mixing of p' with other p-resonances. Suppose one uses a specific dynamical model of 
p resonances predicting the normalization factors c P)P t t .. of the BW-propagators. After in- 
cluding the mixing, the amplitude will have the form (|3*3j) with a complicated s-dependence 
including an imaginary part. This system of mixed propagators could then be diagonal- 
ized, giving in the general case rise to complex couplings between vector mesons and pions. 
Since we do not attempt to solve this complicated dynamical pattern, in the model of our 
choice we will keep the coefficients c P)P i r , for few first resonances as free fit parameters and 
for simplicity, restrict ourselves to real values. 

For the description of the infinite series of higher excitations we adopt an ansatz rooted 
in the Veneziano amplitude and dual resonance models. The specific dual-QCDN c=OD am- 
plitude, suggested in [20], contains an infinite amount of zero-width vector mesons with 
the quantum numbers of p: 

^)=E^. (36) 

For convenience we will count p-resonances by a number n which starts from n = for the p 
meson, so that p'(1450) and p (1700) correspond ton = 1,2, respectively. The coefficients 

(-i) n r(/3 - 1/2) 



a>y/7fm 2 n T(n + l)T([3-l-n) 

decrease rapidly. The parameter a' = l/{2m 2 p ) is related to the p-meson Regge trajectory: 
a p (s) = 1 + a'(s — m 2 p ). Furthermore, the model postulates an equidistant mass spectrum: 

m 2 n = m 2 p (l + 2n) . (38) 

The parameter (3 is free and has to be fitted. We will use Cq fitted from the p region and 
calculate f3 from Eq. (pTTJ). As we shall see, the fit yields c = 1.098 — 1.171 corresponding 
to (3 = 2.16 -T- 2.3, in agreement with [20J. 

Using the above assumptions one easily obtains the form factor in the closed analytical 
form: 

T(B - 1/2) 

FM = r v(R '-k B{p - 1, 1/2 - a's) , (39) 
y/nT((3 - 1) 

where B(x,y) is Euler's Beta-function, so that ^(0) = 1. Importantly, in this model also 
the mean-squared charge radius of the pion (defined as (r£) = 6(dF n (s) / ds)\ s=0 ): 

(rl) = 0.42 4- 0.44 fm 2 , (40) 
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at (3 — 2.16 -T- 2.3, agrees well with the experimental value (301 ( r l)ex P = 0.439 ± 0.008 fm 2 
and with the recent chiral perturbation theory determination |3T] {rfychPT — 0.452±0.013 
fm 2 . 

The most spectacular property of the dual-QCD7v c =oo model (inherited from Veneziano 
amplitude) is its explicit duality: in the timelike region the form factor has poles located 
at s = m 2 , whereas in the spacelike region at large s < it exhibits a smooth behaviour 
with a power-law asymptotic behaviour determined by the parameter (3: 

ton ~ -Z-i • ( 41 ) 



s— >— CO 



With /3 — 2.1 — 2.3 this is very close to the prediction of perturbative QCD. 

In the model was further improved by including the constant widths of resonances 
through the replacement: 



2 2 

mi mi 



(42) 



m 2 — s — ie m\ — s — im n Y n 

The ansatz [2TJ |2U| adopted for the total widths is again motivated by string-like models: 

r„ = 7m„ , (43) 

with 7 = 0.2 fixed from p. As explained above, to account for the presence of 2n and other 
intermediate multiparticle states coupled to p resonances one has to modify the widths 
to include s-dependence with a proper threshold behaviour. Otherwise, the form factor 
predicts an unphysical imaginary part at s = 0. Being unable to account for all possible 
intermediate multiparticle states coupled to each p n we use, as a remedy, the threshold 
behaviour of p n — ► 2n partial width attributing it to the total width 



4 The model for the pion form factor 

After all these modifications the model j2H| for the pion form factor becomes 

oo 

F w { 8 ) = J2 c nBW n (s). (45) 

n=0 

In its simplest form it depends only on few parameters (j3, a' and 7) and includes infinitely 
many hadronic degrees of freedom. 

Importantly, the coefficients c n decrease at n — > 00. Hence starting from n ~ 4, 5, 
moderate deviations of c n ,m n ,T n from the model predictions do not influence the form 
factor, at least in the region of our interest, at ^fs < 2 — 2.5 GeV. On the other hand, 
for the most important first four resonances (including p'" = p% with 777.3 —2 GeV) we 
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will allow the coefficients, masses and widths to deviate from the dual-QCD7v c =oo model 
values, having in mind the effects of coupling to intermediate multiparticle states discussed 
in the previous section. To large extent, the fitted parameters will be close to those of the 
dual-QCDAr c=00 model (see Table HJ). 

Summarizing, our model for the pion form factor has the following form: 



n=0 



+ 

fit 



=4 



(46) 

dual—QCD^ c=oc 



where in the p contribution (n = 0) the p — u mixing is included according to Eq. with 
the fixed parameters taken from [J2ij- In the above, the parameters of the four lowest po,i,2,3 
states (i.e. p, p'(1450), p"(1700) and p'") are supposed to be fitted to experimental data, 
whereas the "tail" with the infinite amount of p n >4 states is taken as in the dual-QCD7v c =oo 
model. However, having in mind the insufficient precision of the current data at y/s > 1 
GeV, we restrict the number of free fit parameters to the coefficients 00,1,2, the masses mo 
and mi and the total widths To and F±. The values of I^ and m2 are taken from 26J. The 
coefficient c 3 of p 3 is fixed from the normalization condition for the form factor: 

c 3 = 1 - (co + ci + c 2 ) fit - I ^2 ° n ) ■ ( 47 ) 

\n=4 / dual-QCD Nc=ao 

All remaining parameters in Eq. (|4^j). that is c n > 4 , m n > 3 , and r n > 3 are calculated from 
Eqs. (jnZ|), (EH) and respectively. 

Let us emphasize the main qualitative features of this model. It nicely matches the 
existing p-dominance models at \fs < 1 GeV, such as the ones considered in ^2], simply 
because the p n >2 states play a minor role in the \/s < 1 GeV region. The model ()46|) is 
flexible, that is, it allows to vary the proportion of fitted and modelled resonances above 
p. E.g., with sufficiently precise data at higher energies one can include also the p4-state 
into the 'fit' part. Alternatively, p 3 can be removed from the 'fit' part and added to the 
dual-QCDTv^oo part. Furthermore, as mentioned already, F n (s), as given by Eq. (|4b|). 
can be easily continued to s < and compared with the experimental data and QCD 
predictions in the spacelike region. Accordingly, one gets a smooth power-like behaviour 
at asymptotically large timelike s. Altogether, the model contains a reasonably small 
amount of free parameters: three per each fitted resonance (the mass, coefficient and total 
width) and three "global" parameters 7, with a' = l/{2m 2 p ) taken from the Regge 
trajectory, j3 taken from Eq. (J37j) with cq derived from the fit, and with 7 = 0.2. In 
principle, one can try to relax and independently fit also the three "global" parameters, 
but we prefer to keep them fixed. Needless to say, the suggested ansatz has considerable 
room for improvement, especially concerning the treatment of the effective s-dependent 
resonance widths. 

We have fitted the model (J4T)|) to the existing data El EE El for the timelike pion form 
factor. The results of the fit for the two different versions of BW-formulae: KS (Eq. (J3U|) 
for all resonances) and GS (Eq. ()31|) for the first three resonances po,i,2), together with the 
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relevant input parameters are presented in Table ^ and compared with the predictions of 
the dual-QCD7v c=0 o model. The results for (F^s)! 2 are plotted in Fig. 1, separately for the 



Pq t*q tti pf pr 


Tnnn t 


FitfKS") 


FitCGSl 

J- 1 U I V.J LJ J 


mi p\ 1_ 


PDG value 










U(jU /V — tYi 




p 




773.9 ±0.6 


776.3 ±0.6 


input 


775.5 ±0.5 


p 




144.9 ±1.0 


150.5 ±1.0 


input 


150.3 ±1.6 




783.0 








782.59 ±0.11 




8.4 








8.49 ±0.08 


Uln' 




1357 ±18 


1380 ±18 


1335 


1465 ± 25 


p 




437 ±60 


340 ± 53 


266 


400 ± 60 


rripii 


1700 






1724 


1720 ±20 


Tpn 


240 


- 


- 


344 


250 ±100 


Trip'" 








2040 




1 ,-.ttt 
p 








400 




c 




1.171±0.007 


1.098±0.005 


1.171 




p 


c and Eq. (jSTJ) 


2.30±0.01 


2.16±0.015 


2.3(input) 






0.00184(KS) 












0.00195(GS) 










Cl 




-0.119±0.011 


-0.069 ±0.009 


-0.1171 




c 2 




0.0115 ±0.0064 


0.0216 ±0.0064 


-0.0246 




c 3 


Eq. © 


-0.0438^0.02 


-0.0309^0.02 


-0.00995 




oo 

71=4 


-0.01936 






-0.01936 




rVd.o./. 




155/101 


153/101 







Table 1: Parameters of the pion form factor \46] ) and results of the fit to the data. Masses 
and widths are given in MeV. The row 'Fit KS( GS) ' contains the fitted values for the case 

oo 

where all BW n are taken as in Eq. A3U\) (BWqi2 taken as in Eq. The sum E c n 

71=4 

is calculated from Eq. \3T\) . The PDG parameters for p(770) are those listed in for 
"Charged only, r decays and e + e~ ". The parameter c w is taken from 11^ - 

\fs < 1 GeV and \fs > 1 GeV regions. 
A few comments are in order. 

The most spectacular result of the fit is the change of the p" coefficient c 2 with respect 
to the model (201 anc ^ to the earlier fits ^2]; from negative values ~ —(0.02 — 0.04) to 
smaller but positive values ~ 0.01-0.02. The positive sign is a direct consequence of the 
dip in the cross section around 1.6 GeV (also the earlier fit of the data in [7] revealed a 
similar pattern). 

Furthermore, the fitted mass of p' gets shifted with respect to the PDG value, the 
latter obtained by adding together data from all decay channels of p' . Note that a lower 
rripi consistent with our fit is also obtained by the fits in [TJ] and predicted by the dual- 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v^[GeV] 

(a) 



10 r 




0.0001 1 1 — 1 u 1 1 1 

1 1.5 2 2.5 3 3.5 

v^[GeV] 

(b) 

Figure 1: The pion form factor squared |F„-(s)| 2 as a function of y/s fitted to the data in 
the region near (a) and above (b) p resonance. The solid (dashed) line corresponds to the 
KS(GS) parameterization of BW formula. Data are taken from 0, '(crosses),^ 1 (stars), 
and (squares) . The triangle point is the value of the form factor extracted from J/ip — ► 
2tt. 
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QCDat c=00 model. It is quite probable that a more elaborated model of the total width of 
p' including multiparticle thresholds will increase the fitted mass. We note that the values 
of the masses and widths of p, p', p" as well as cq and c\ are in the ballpark of the dual- 
QCDat c=00 model. The same is true for the magnitude of C2, its positive sign is enforced 
by the dip around 1.6 GeV and might be a consequence of the strong mixing between p" 
and nearby resonances. 

In general, the p', p" and p" terms and their interplay with the contribution of p in 
both imaginary and real parts of the form factor are the main effects which determine the 
behavior of |F„-(s)| 2 at yfs > 1 GeV. In particular, the dip in the form factor observed 
near y/s = 1.6 GeV can only be described by altering the sign of c^. The role of the 
summed "tail" of p n >\ states is less important. One has to admit that the quality of the 
fit is not very high, (we get typically x 2 /d-o.f. ~ 1.5), in fact this could simply indicate 
inconsistencies in the normalization of the various pieces of existing data in the timelike 
region. Moreover, the data points with large errors at y/s = 2.5 4- 3 GeV are systematically 
higher than the fitted curve, and the value of F 7r (y/s = mj/^) extracted from the J/ip — > 2n 
partial width is larger than the model prediction by a factor of about three (this point was 
not included into the fit). This disagreement deserves a comment. Note that the extraction 
of F v (m 2 j,^) is "theoretically biased", because one tacitly assumes that the intermediate 
photon exchange is the only mechanism in this decay. Other mechanisms such as one- 
photon plus two intermediate gluons could also be important in the J /if} — > 2n. Hence, the 
hadronic matrix element in this isospin-violating transition could be actually different from 
the pion form factor. Although estimates [29| of gluonic effects based on the perturbative 
charmonium annihilation are in favour of their smallness, we still think there is room for 
nonperturbative effects which are however not easily assessed. For the time being it seems 
difficult to accommodate a pion form factor as large as the one derived from J/if> decay. 
(It will be interesting to check for a similar effect in the 47r mode.) 

To check our model further we continue F n (s) to s < and compare the result with 
the data there (see Fig. 2). The form factor obtained from direct electron-pion scattering 
at small s < (up to \s\ = 0.25 GeV 2 ) |3U] is not sensitive to the contributions of higher 
than p resonances, provided the normalization to one at s = is imposed. Note that the 
fitted form factor (J4T)J) predicts (r 7r ) 2 = 0.440(0.426) fm 2 for the KS(GS) versions, close to 
the dual-QCD7v c=0 o prediction pUj) and to the experimental value jSOj quoted above. 

At larger \s\ the old data [32] obtained from pion electroproduction have large errors 
and suffer from some intrinsic uncertainties (SHI- More accurate data obtained recently 
at JLab [31] at \s\ < 1.6 GeV are in agreement with the model, but not sensitive to the 
details of the fit. Furthermore, there is reasonable agreement between the model and the 
QCD light-cone sum rule (LCSR) predictions The latter are taken from ^7] with 

their estimated theoretical uncertainty. The pion form factor calculated from 3-point QCD 
sum rules ^3] in the region of their validity \s\ ~ 1 — 4 GeV 2 is within the LCSR interval 
and therefore not shown separately. 

If we try to artificially enhance the form factor at large timelike region, e.g., by en- 
hancing the contribution of the "tail" trying to fit also the point at yfs = mj/^, the form 
factor at spacelike s < increases correspondingly, and the general agreement between 
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(b) 

Figure 2: The pion form factor squared \F n (s)\ 2 as a function of s in the spacelike region at 
low \s\ (a) and large \ s\ (b). The solid (dashed) line corresponds to the analytic continuation 
of the timelike form factor with the KS(GS) parameterizations of the BW formula. 
Data are taken from pity (crosses). pDj^ ( open circles) and pfltf (full squares). Dotted lines 
at \s\ > lGeV 2 represent the interval derived from QCD light-cone sum rule predictions 



data and QCD sum rule predictions is lost. We conclude therefore, just as before, that 
it is implausible for the form factor obtained on the basis of dual resonance approach to 
reach values l-F^s)) 2 > 0.01 at <Js = 2.5 4- 3 GeV without conflicting with the spacelike 
data and especially with QCD predictions. This statement is independent of many details 
involved in the timelike form factor model and in the QCD calculations. It is therefore 
extremely interesting to obtain new accurate data at least up to \fs = 2.5 GeV to check 
this conjecture. 



5 Charged and neutral kaon e.m. form factors 

We now adopt the analogous strategy to describe the kaon form factors. Combining in- 
formation on K + K~ and K°K° production with constraints from isospin symmetry, and 
using assumptions deduced from the quark model and the OZI rule, it is possible to sepa- 
rate the 1 = 1 and 1 = amplitudes in the form factor, and even the u- and ^-components 
of the 1 = part. The 1=1 part can then be used to predict the rate for r — ■> vK~K°. 

The electromagnetic form factors for charged and neutral kaons defined similar to 
Eq. (P): 

(K+( Pl )K-(p 2 )\j™\0) = ( Pl -p 2 ),F K+ (s) (48) 
(K°( Pl )K (p 2 )\jjr\0) = (p 1 -p 2 ),F K o(s) , (49) 

obey the constraints 

F K+ (0) = 1, F K o(0)=0. (50) 
They can be separated into their isospin one and zero parts, 

J"K+(K°) ~ ^ K +(X0) + K+(K )' l 5i ) 

From isospin invariance one derives 

and the 1 = 1 part can furthermore be used to predict the charged current matrix element: 

(K + (Pi)K°(p 2 )\j,\0} = ( Pl -p 2 )»2F«f\s) . (53) 

A simultaneous fit to the two electromagnetic form factors leads, therefore, to a direct 
prediction for the rate for r — > vK~K°, without any further assumption. 

In the context of vector dominance, combined with the quark model, the kaon form 
factors are saturated by p, uj, (ft and their radial excitations, 

F K (s) = > — ; — , 54) 

v Ti , a, Wy-s- im v Y v 

V=p,u),q>,p' ,4> v 

and it is the p-mediated 1=1 part which enters both the e.m. and the charged current 
matrix elements. We define the decay constants of the vector mesons via 

(V|j7|0) = Kvmvfvfa ' ( 55 ) 
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A 



qs 



A 



sq 



Figure 3: The quark diagrams corresponding to the contributions to the pion and kaon form 
factors with various flavour combinations. 

where ey is the polarization vector of V, and the coefficients k p = l/y/2 (see Eq. (0)), k w = 
l/(3\/2) and = —1/3 reflect the valence quark content of these mesons corresponding 
to "ideal" mixing: 



P 



uu — dd 



uu + dd 



10 



ss . 



V2 ' ~ y/2 

The strong coupling is defined as in Eq. (jlOjl : 

(K( Pl )K(p 2 )\V) = { P 2-pi) v e v[v) g VKR . 
In the flavour SU(3)-symmetry limit one evidently has 

F K+ (s)=F 7V (s), F K o(s) = 0. 



(56) 



(57) 



(58) 



Since this symmetry is broken by the strange-nonstrange quark-mass difference, one has to 
expect quite noticeable deviations from Eq. (|58j) . in particular the K° form factor does not 
vanish at nonzero s. Moreover, since BR(<fi — > K + K~) ~ BR{<p — > K°K°), the neutral 
and charged kaon form factors have almost equal magnitudes near the <ft resonance. In 
fact, the SU(3)-breaking also manifests itself in the valence quark content of vector mesons 
given in Eq. (|5T)|) and in the mass splitting between m p ~ and m^. 

In order to derive the kaon form factors in terms of separate vector meson contribu- 
tions, it is convenient to consider a generic quark diagram of the strong VPP coupling 
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and distinguish these diagrams by the presence and position of s quarks. In the isospin 
symmetry limit, which we adopt here, there are three different diagrams depicted in Fig. 3: 
1) without strange quarks (upper), 2) with s and s in the PP state only (middle) and 3) 
with s and s in the V and the PP state (bottom). We denote the corresponding hadronic 
invariant amplitudes A qq , A qs and A sq . The diagrams with charge conjugated quark lines 
lead to the same amplitudes with an additional minus sign, taking into account the nega- 
tive C-parity of neutral vector mesons. The strong couplings of p, u, are then expressed 
in terms of diagrams: 

9p°Tr + w~ = 9pmr V^2^4gq j 

1 A 1 A 

9pOK+K 9uK+K — ^= A qs , 9 p K 0K0 — —g^^RO ~ — A qs , 

9(pK+K- — 9<t,K°K = ~A sq . (59) 

Thus, from the quark model one expects approximately equal KKp and KKuj couplings. 
It is easy to check that this simple formalism correctly reproduces SW+tt- = g^o^o = 0, 
as well as SU(3)-symmetry relations between the couplings. In what follows we also use 
the relation between the decay constants of u and p following from Eq. (jB3|) : j w = f p . 
Substituting the decay constants and hadronic couplings (JoT^) to Eq. (p>^) we obtain the 
desired decompositions of the kaon form factors in terms of vector meson contributions: 



F,As) = -k^BW^s) + i^BWUs) + ! -i^BW t ( S ) . (61) 
Imp brn^ Srritf, 

Written in the same terms the pion form factor is 

Fn (s) = ^BW p (s) , (62) 
m p 

and f p A qq /mp = 1 in the simplest version of VDM. In the SU(3) limit f p = f^, m p = = 
, A qq = A qs = A sq and the relations (Jo"H|) are reproduced. In Eqs. (JBTiJ) and (joT|) we will 
adopt the KS version of BW formulae for p, and the analogous s-dependent width for 0, 



ml ( p K (s) 



p K {ml) 



3 

r . (63) 



with = r to *(0) and p K (s) = (s — Am 2 K ) l l 2 /2. For simplicity, we assume that the effective 
threshold of all decay modes including — > 3tt is approximated by Eq. ()63|). having in 
mind that the non- KK channels give only about 20% of T^. Naturally, the possibility to 
use the GS-form in Eqs. (160)) and (|61|) exists, yielding inessential differences. For simplicity, 
to avoid complicated 37r-threshold factors we will use constant widths for u, having in mind 
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that the thresholds are much lower than the boundary of the physical region of the form 
factor: 9m^ <C Am 2 K . This is a good approximation at least for the narrow uo resonance. 

Adding radial excitations to all ground-state vector mesons is the next natural step. 
From the pion form factor analysis we already learned that the "tails" of higher resonances 
are numerically inessential. For the kaon form factor we therefore restrict the analysis to 
the excited states p', p", J = cj(1420) , lu" = cj(1650) and 0' = 0(1680) Higher 
excitations, as well as more elaborated s-dependent widths, can be installed in the future 
when more accurate data will be available. Since the products of decay constants and 
strong couplings in Eqs. (f6T)|) and (f6T]) will not be separated and have to be fitted as a 
whole, it is convenient to introduce again the normalization factors Cy instead of these 
products. The ansatz for the kaon form factors thus reads: 



Fk + (s) = ^BW p {s) + c«BW p ,{s) + $,BW p „{s)) 
1 
6 



+l(c£BW u (s) + cfiBW^s) + C k,BWu»(s)) 



+i(c*SW*(s) + c^BW^(s)) , (64) 

F K o{s) = - l -(c K p BW p (s) + cf,BW p ,(s) + cf„BW p „(s)) 

+hc*BW u {s) + c^BW^s) + c^BW^is)) 
6 

+J( W Bl^(s) + CffBW^s)) , (65) 

The widths are with p-wave factors for p and (f) states as explained above, and constant for 
a»-states, which is however a rather crude approximation for lu', lu". The ansatz in Eqs. (JMj) 
and (J65|) reflects isospin invariance and the hierarchy of vector meson contributions accord- 
ing to their valence-quark content, however, it allows for the possibility of SU(3) violations 
which could and will manifest in differences between the fitted normalization coefficients. 
The additional factor in Eq. (J55J) takes into account the isospin-breaking difference 
between the charged and neutral kaon couplings to 0: 

= g.KORQ = / ^(0^^°)K-4m^ + p \ 1/2 
W " g*K +K - \BR(<j> - K+K-)(m* - J " TO 

According to the central value of this factor slightly deviates from the unit: 

^ = 1.027 ±0.01. (67) 

In the vicinity of the <fi resonance this small effect is noticeable in the fit, and as far as the 
branching ratio is concerned, is dominated by the phase space factor. The factor rj^ also 
takes care of Coulomb-rescattering and other isospin-violating differences between charged 
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and neutral modes (see also [HSHHS])- To ensure the proper normalizations F K o(0) = and 
FjK+ify = 1> we introduce an additional energy dependence with a simple step-function 

V<j> {s) = 1 + (ty - - (m - T^eim^ + T (t> -^), (68) 

which in future, after this effect is better understood both experimentally and theoretically, 
can be replaced by an appropriate analytical energy-dependence. 

We have fitted the model (JH1 and §B> to the available data on charged jS3 EH1 EH1 113 
l4~T] and neutral [H EH El SB] kaon form factors. The masses and widths of p, uj and their 
excitations are taken from (20] and are listed in the Table EJ Two different variants of the 
fit are carried out: 

(1) the constrained fit (motivated by the quark model) where the normalization factors 
for uj resonances are fixed: c% u , w „ = cf p ' <p " and only the normalization factors for the p 
resonances are fitted; 

(2) the unconstrained fit, where uj- and p- factors are fitted as independent parameters. 
First, the mass and width of <fi as well as the coefficient r\§ are fitted in the region 

around resonance. We obtain: r]^ = 1.011 ± 0.009 (1.019 ± 0.009 ) for the constrained 
(unconstrained) fit, in a good agreement with the experimental value JB7J). Fixing r]^ 
and using and as starting values, the data in the whole region of ^/s are then 
fitted. The results of the fit are collected in Tabled The best (i.e., stable and physically 
plausible) results for both variants of the fit are obtained if data on and F K o are 
fitted simultaneously. Thus, predicting F K o from F K + with the currently available data is 
not yet possible. The resulting curves for the form factors are plotted in Figs. 4,5. Most 
importantly, fitting the kaon form factor above resonance, it is indeed possible to extract 
separate p,uj,<p components, which was not possible in the <fi region due to the dominance 
of this resonance. 

We also find the pattern of the normalization factors c^ p , for the first two p-resonance 
to be very similar to the corresponding values Co t i obtained in the pion form factor fits. 
These factors can be immediately translated into the strong couplings dividing out the 
decay constants of vector mesons. The latter are independently measured in the leptonic 
decays revealing a very mild SU(3) breaking, at the level of 10 %; according to the data 
in |2E]: f p ~ 220 MeV, f u ~ 195 MeV and / ~ 228 MeV. The SU(3)-violating difference 
between the couplings of p and (p to kaons estimated from comparing and c^ is also 
moderate in both versions of the fit. 

As already noticed above, the constraint c^ = c^ naturally follows from the valence 
quark content of both mesons and we consider this constraint as a part of our model. The 
fact that the unconstrained fit gives about 25% difference between these two coefficients, 
a noticeable deviation from the quark-diagram relation, should be taken with caution, 
having in mind poor quality of data. Also x 2 ' s of both fits are in the same ballpark, 
so that from the fitting point of view we cannot yet give any preference to the version 
with the "floating" couplings of u;-resonances. On the other hand, this difference indicates 
that the fit is able to resolve also the "fine structure" of the couplings. The differences 
between the normalization factors given by the fit for excited resonances: Cy vs. c^, (in the 
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Parameter 


Input 


Fit(l) 


Fit(2) 


PDG value 26 




_ 


1019.372 ±0.02 


1019.355 ±0.02 


1019.456 ±0.02 




_ 


4.36 ±0.05 


4.29 ±0.05 


4.26 ±0.05 




1680 


_ 


_ 


1680 ±20 




150 


_ 


_ 


150 ±50 


Vtlp 


775 






775.8 ±0.5 


r p 


150 


_ 


_ 


150.3 ±1.6 




1465 


- 


- 


1465 ± 25 




400 


- 


- 


400 ± 60 


77lp// 


1720 


_ 


_ 


1720 ±20 


Y p ii 


250 


- 


- 


250 ±100 


m w 


783.0 


- 


- 


782.59 ±0.11 


K 


8.4 


- 


- 


8.49 ±0.08 




1425 


- 


- 


1400-1450 


rv 


215 


- 


- 


180-250 




1670 


- 


- 


1670 ± 30 




315 


- 


- 


315 ±35 


C0 




1.018 ±0.006 


0.999 ±0.007 




C0' 




-0.018 =F 0.006 


0.001^0.007 








1.195 ±0.009 


1.139±0.010 








-0.112±0.010 


-0.124±0.012 








-0.083^0.019 


-0.015 T 0.022 






c p 


1.195 ±0.009 












1.467 ±0.035 




<£(!) 


c p' 


-0.112±0.010 






<#(2) 






-0.018 ±0.024 






1 _ r K _ r K 

1 C w C UJ' 


-0.083^0.019 


-0.449^0.059 








328/242 


281/240 





Table 2: Parameters of the kaon form factors and results of the fit to the data. Masses and 
widths are given in MeV. The row 'Fit(l) 7 (Fit(2)) contains the values of the constrained 
(unconstrained) fits. 
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unconstrained fit), c5 vs. c^, etc. are generally large, which is not surprising, in view of 
the complicated mixing between all these states. Including in the future more precise data 
and switching on the "tails" of the dual QCDat c=00 amplitudes in the kaon form factors for 
all three vector mesons will allow to reveal these differences more accurately. 

Furthermore, an indication for an excess of the measured charged kaon form factor vs. 
the model is present in Fig. 4b in the region around 2 GeV, although the experimental 
errors are large. Remember, that we have not included in our fit the contribution of the 
second excited </>" state with a mass around 2 GeV, which might be responsible for this 
potential difference. However, we refrain from further investigation before more accurate 
data are available. 

The mean-squared charge radius of the K + obtained in our model: a/ (r 2 K ) = 0.56 fm 
(for both fits and with a small error), is in a good agreement with the experimental value 
jSj \J ( r K) e xp = 0.53 ±0.05 fm. We have also checked that, being analytically continued to 
large s < — lGeV 2 , the charged kaon form factor agrees with the LCSR prediction obtained 
in H2|. 

Finally, as mentioned in the Introduction, the separate reconstruction of / = and 
1=1 spectral functions might be a useful ingredient for various phenomenological analyses. 
In Fig. 6 we display the spectral functions defined as 



12tt 



±F K o(s) 



2p K (s) 



(69) 



noticing that this observable is quite sensitive to the pattern of resonances in the form 
factor. 



6 Predicting r — > K K°u T decay distribution and rate 

As emphasized above, the isospin one part of the e.m. kaon form factor, together with the 
isospin-symmetry relation (J53j) . can be used to predict the r — > K~K°i/ T decay width. The 
differential decay distribution in \[Q 2 (the invariant mass of the kaon pair), normalized to 
the leptonic width of r reads: 

1 ^ dBR(r -> K~K°v T ) 

)2 



BR(t -> yrvyVr) ) d^J~Q 



2ml \ ml J \ ml J \ Q 

In accordance with the isospin limit, we neglect the mass difference between charged and 
neutral kaons and the effect of the scalar form factor. Using Eq. (J53j) we have F K - K o = 
-2F^+ 1] , hence 

\F k -k°(Q 2 )\ 2 = \c«BW p {Q 2 ) + c K p ,BW p ,{Q 2 ) + c«BW p „(Q 2 )\ 2 . (71) 
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Figure 4: The charged kaon form factor squared \F K +(s)\ 2 as a function of y/s fitted 
to the data taken from fff?| / (crosses), pfitf (open squares), fffffi (open circles), J%U/ (full 
squares) and 141)1 (full circles). The solid(dashed) lines correspond to the constrained (un- 
constrained) fit. 
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Figure 5: The neutral kaon form factor squared \Fxo(s)\ 2 as a function of \/s fitted to the 
data taken from ^(triangles), \3l\l (crosses), \4M ( open squares) and \J$ (full circles). The 
solid(dashed) lines correspond to the constrained (unconstrained) fit. 
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Figure 6: The spectral functions \69) with 1 = (a) and 1=1 (b) obtained from the fitted 
kaon form factors. The solid(dashed) lines correspond to the constrained (unconstrained) 
fit. 
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Figure 7: The normaized distribution t{t^k-k {) u ) — ^ n ^ e ^ aon V a ^ r invariant mass 

\[Q^ in units of GeV' 1 obtained from the fitted kaon form factor; the solid (dashed) line 
corresponds to the constrained (unconstrained) fit. The event distribution measured by 
CLEO Collaboration [45Jj and normalized, dividing by the total number of events, is shown 
with points. 



The fitted values for c^ pl „ from Table EJ thus allow us to calculate the decay distribution 
(I7(J|) . The normalized distribution is plotted in Fig. 7 and (qualitatively) compared with the 
event distribution in the kaon pair mass, measured by CLEO Collaboration |45j . Another 
measurement of this distribution by ALEPH Collaboration can be found in ^B] • Integrating 
over \[Q^ from 2m^ to m T we obtain the branching ratio 

BR(t -> K-K°u T ) = 0.19 ±0.01% (0.13 ±0.01%) (72) 

for the constrained (unconstrained) fit, to be compared with the experimentally measured 
value pp] 

BR(t^K-K°v t ) =0.154 ±0.016%. (73) 

We see that both the decay distribution and the decay width are very sensitive to the 
pattern of p resonances in the isospin-1 form factor. Generally, the width grows with the 
increase of the excited p contributions, an effect observed earlier in : 4 .7] (see also jH]). 
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7 Conclusions 



In this paper we considered the models of timelike form factors of pions and kaons in 
anticipation of new and more accurate data in the region above 1 GeV, from e + e~ ma- 
chines using the radiative return method. We introduced an ansatz for the pion form 
factor which is based on dual-resonance models and Veneziano amplitude. We argued that 
the parameters of the ground-state and first excited states can deviate from the model 
prediction due to effects of mixing with multiparticle (e.g., two - or 4 tt states), therefore 
have to be fitted independently as free parameters. From the fit to the available pion form 
factor data we have found that the main contribution to the form factor originates from 
the ground state plus the first two-three radially excited states. The tail from the infinite 
series of resonances produces inessential, but visible effects. The sign and value of the coef- 
ficients at certain excited resonances are shifted with respect to the dual QCDjv c=00 ansatz 
signaling large mixing effects. Possible checks of the model are provided by the spacelike 
form factors, the pion charge radius and the behaviour for large s in the timelike region. 
In particular around y/s = 3 GeV we predict a value smaller than the one anticipated from 
J I ip decay. On the other hand, data fitted to this model can be used for important tests of 
quark-hadron duality, and of the QCD calculations of the pion form factor in the spacelike 
region. 

Furthermore, we formulated an analogous model for the kaon form factor and demon- 
strated that the contributions of <fi , u and p resonances (or, alternatively, the isospin 
and 1 components) can be separated by the fit. Interestingly, the r — > K~ K°i> T -deceLy 
distribution and partial width predicted from the model manifest a substantial sensitivity 
to the pattern of p resonances in the isospin-1 part of the form factor. 

The model still has considerable room for improvement. In particular, a more detailed 
kinematical and dynamical analysis of total widths in the Breit-Wigner factors would allow 
to implement a more accurate energy-dependence in these widths. 
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